For the past few decades, the scientific community has been locked in a fierce debate over the issue of climate change. Those concerned about the future of the environment, known as “alarmists,” warn that increased carbon dioxide emissions and the resulting global temperature increase of about 1.7 degrees Fahrenheit since 1880 (NASA) will have disastrous effects on agriculture, ecosystems, water resources, and human health. With increased temperatures come the risks of heat waves, drought, early snow melt, and extreme weather events (EPA). Although “skeptics” may acknowledge increased temperatures, they deny that there is anything out of the ordinary that should cause humans to change behavior.

The Golden State of California, already known for its sunny skies and breezy coasts, has been warming for the past century. Southern California specifically has experienced a dramatic temperature increase of about 3 degrees Fahrenheit in the past century. If this is the case, then why have the maximum temperatures (TMax) of more than half the year in Ojai, California decreased since 1920?

Ojai is a small town of about 7500 people situated in Ventura County, just northwest of Los Angeles and east of Santa Barbara. According to National Oceanic and Atmospheric Association climate data, the maximum temperatures in Ojai have decreased for seven out of twelve months of the year. Out of these seven months, five have p-values less than 0.05. Therefore this data is statistically significant, which means that it rejects the null hypothesis that modern climate change is due to natural causes. Four of the statistically significant months have decreasing temperatures.

## [1] "19170101" "19170102" "19170103" "19170104" "19170105" "19170106"
## [1] OJAI CA US
## Levels: OJAI CA US
## 
## Call:
## lm(formula = TMAX ~ NewDate, data = LosAngeles)
## 
## Residuals:
##     Min      1Q  Median      3Q     Max 
## -44.919  -9.415  -0.342   9.350  40.416 
## 
## Coefficients:
##               Estimate Std. Error  t value Pr(>|t|)    
## (Intercept)  7.807e+01  6.640e-02 1175.888  < 2e-16 ***
## NewDate     -2.658e-05  6.259e-06   -4.247 2.17e-05 ***
## ---
## Signif. codes:  0 '***' 0.001 '**' 0.01 '*' 0.05 '.' 0.1 ' ' 1
## 
## Residual standard error: 12.41 on 35529 degrees of freedom
##   (782 observations deleted due to missingness)
## Multiple R-squared:  0.0005074,  Adjusted R-squared:  0.0004793 
## F-statistic: 18.04 on 1 and 35529 DF,  p-value: 2.172e-05

##   Month Year      TMAX YEAR MONTH  NewDate
## 1    01 1917  7.069547 1917     1 1917.000
## 2    02 1917  8.784908 1917     2 1917.083
## 3    03 1917  7.809974 1917     3 1917.167
## 4    04 1917  8.435284 1917     4 1917.250
## 5    05 1917  6.535758 1917     5 1917.333
## 6    06 1917 12.658194 1917     6 1917.417

## 
## Call:
## lm(formula = TMAX ~ YEAR, data = MonthlyMean[MonthlyMean$Month == 
##     "04", ])
## 
## Residuals:
##      Min       1Q   Median       3Q      Max 
## -10.1433  -2.5047  -0.0938   2.8245   8.0180 
## 
## Coefficients:
##             Estimate Std. Error t value Pr(>|t|)  
## (Intercept)  21.6564    25.9442   0.835   0.4059  
## YEAR          0.0267     0.0132   2.023   0.0458 *
## ---
## Signif. codes:  0 '***' 0.001 '**' 0.01 '*' 0.05 '.' 0.1 ' ' 1
## 
## Residual standard error: 3.75 on 96 degrees of freedom
## Multiple R-squared:  0.0409, Adjusted R-squared:  0.03091 
## F-statistic: 4.094 on 1 and 96 DF,  p-value: 0.04582

## 
## Call:
## lm(formula = TMAX ~ YEAR, data = MonthlyMean[MonthlyMean$Month == 
##     "07", ])
## 
## Residuals:
##     Min      1Q  Median      3Q     Max 
## -8.9677 -2.0630 -0.1102  2.3523  8.7120 
## 
## Coefficients:
##              Estimate Std. Error t value Pr(>|t|)    
## (Intercept) 191.52739   22.20627   8.625 1.24e-13 ***
## YEAR         -0.05114    0.01129  -4.529 1.69e-05 ***
## ---
## Signif. codes:  0 '***' 0.001 '**' 0.01 '*' 0.05 '.' 0.1 ' ' 1
## 
## Residual standard error: 3.238 on 97 degrees of freedom
## Multiple R-squared:  0.1745, Adjusted R-squared:  0.166 
## F-statistic: 20.51 on 1 and 97 DF,  p-value: 1.69e-05

## 
## Call:
## lm(formula = TMAX ~ YEAR, data = MonthlyMean[MonthlyMean$Month == 
##     "08", ])
## 
## Residuals:
##     Min      1Q  Median      3Q     Max 
## -6.8232 -1.5971 -0.0065  1.5255  8.5973 
## 
## Coefficients:
##              Estimate Std. Error t value Pr(>|t|)    
## (Intercept) 134.19879   20.15496   6.658 1.61e-09 ***
## YEAR         -0.02162    0.01025  -2.109   0.0375 *  
## ---
## Signif. codes:  0 '***' 0.001 '**' 0.01 '*' 0.05 '.' 0.1 ' ' 1
## 
## Residual standard error: 2.958 on 98 degrees of freedom
## Multiple R-squared:  0.04343,    Adjusted R-squared:  0.03367 
## F-statistic: 4.449 on 1 and 98 DF,  p-value: 0.03746

## 
## Call:
## lm(formula = TMAX ~ YEAR, data = MonthlyMean[MonthlyMean$Month == 
##     "11", ])
## 
## Residuals:
##     Min      1Q  Median      3Q     Max 
## -6.8244 -3.3259  0.3259  2.7188 10.2763 
## 
## Coefficients:
##              Estimate Std. Error t value Pr(>|t|)    
## (Intercept) 163.69455   27.98160    5.85 6.65e-08 ***
## YEAR         -0.04539    0.01423   -3.19  0.00192 ** 
## ---
## Signif. codes:  0 '***' 0.001 '**' 0.01 '*' 0.05 '.' 0.1 ' ' 1
## 
## Residual standard error: 4.08 on 97 degrees of freedom
## Multiple R-squared:  0.09493,    Adjusted R-squared:  0.0856 
## F-statistic: 10.17 on 1 and 97 DF,  p-value: 0.001919

## 
## Call:
## lm(formula = TMAX ~ YEAR, data = MonthlyMean[MonthlyMean$Month == 
##     "12", ])
## 
## Residuals:
##     Min      1Q  Median      3Q     Max 
## -7.9694 -3.1011 -0.3825  3.0747  9.2401 
## 
## Coefficients:
##              Estimate Std. Error t value Pr(>|t|)    
## (Intercept) 137.85146   27.48570   5.015 2.42e-06 ***
## YEAR         -0.03572    0.01398  -2.555   0.0122 *  
## ---
## Signif. codes:  0 '***' 0.001 '**' 0.01 '*' 0.05 '.' 0.1 ' ' 1
## 
## Residual standard error: 3.958 on 96 degrees of freedom
## Multiple R-squared:  0.06365,    Adjusted R-squared:  0.0539 
## F-statistic: 6.526 on 1 and 96 DF,  p-value: 0.0122

##   Month Year     TMAX YEAR
## 1    01 1917 61.61290 1917
## 2    02 1917 66.71429 1917
## 3    03 1917 69.06452 1917
## 4    04 1917 71.53333 1917
## 5    05 1917 71.87097 1917
## 6    06 1917 90.33333 1917

Interestingly, while Ojai is experiencing a cooling trend in TMax for a majority of months of the year, there is an opposite phenomenon occurring with minimum temperatures (TMin). For 11 out of 12 months, Ojai’s minimum temperatures were increasing. 8 out of these 11 months were statistically significant, with p-values less than 0.05. The one month, December, with decreasing TMin was not statistically significant.

TMIN

##   Month Year     TMIN YEAR
## 1    01 1917 34.58065 1917
## 2    02 1917 35.92857 1917
## 3    03 1917 33.35484 1917
## 4    04 1917 39.46667 1917
## 5    05 1917 41.09677 1917
## 6    06 1917 51.60000 1917

## 
## Call:
## lm(formula = TMIN ~ YEAR, data = MonthlyMeanTMIN[MonthlyMeanTMIN$Month == 
##     "01", ])
## 
## Residuals:
##     Min      1Q  Median      3Q     Max 
## -8.5486 -1.7755 -0.1755  1.3592  7.8665 
## 
## Coefficients:
##              Estimate Std. Error t value Pr(>|t|)   
## (Intercept) -26.78993   19.92827  -1.344  0.18195   
## YEAR          0.03183    0.01013   3.142  0.00222 **
## ---
## Signif. codes:  0 '***' 0.001 '**' 0.01 '*' 0.05 '.' 0.1 ' ' 1
## 
## Residual standard error: 2.948 on 98 degrees of freedom
## Multiple R-squared:  0.09151,    Adjusted R-squared:  0.08224 
## F-statistic: 9.871 on 1 and 98 DF,  p-value: 0.00222

## 
## Call:
## lm(formula = TMIN ~ YEAR, data = MonthlyMeanTMIN[MonthlyMeanTMIN$Month == 
##     "03", ])
## 
## Residuals:
##    Min     1Q Median     3Q    Max 
## -5.073 -1.373 -0.519  1.671  6.960 
## 
## Coefficients:
##               Estimate Std. Error t value Pr(>|t|)    
## (Intercept) -21.363213  17.456061  -1.224 0.224012    
## YEAR          0.031190   0.008879   3.513 0.000678 ***
## ---
## Signif. codes:  0 '***' 0.001 '**' 0.01 '*' 0.05 '.' 0.1 ' ' 1
## 
## Residual standard error: 2.523 on 96 degrees of freedom
## Multiple R-squared:  0.1139, Adjusted R-squared:  0.1047 
## F-statistic: 12.34 on 1 and 96 DF,  p-value: 0.000678

## 
## Call:
## lm(formula = TMIN ~ YEAR, data = MonthlyMeanTMIN[MonthlyMeanTMIN$Month == 
##     "05", ])
## 
## Residuals:
##     Min      1Q  Median      3Q     Max 
## -5.6486 -1.6267  0.0623  1.6795  6.4054 
## 
## Coefficients:
##              Estimate Std. Error t value Pr(>|t|)   
## (Intercept) -5.887086  16.532669  -0.356  0.72256   
## YEAR         0.027011   0.008409   3.212  0.00179 **
## ---
## Signif. codes:  0 '***' 0.001 '**' 0.01 '*' 0.05 '.' 0.1 ' ' 1
## 
## Residual standard error: 2.39 on 96 degrees of freedom
## Multiple R-squared:  0.09704,    Adjusted R-squared:  0.08764 
## F-statistic: 10.32 on 1 and 96 DF,  p-value: 0.001794

## 
## Call:
## lm(formula = TMIN ~ YEAR, data = MonthlyMeanTMIN[MonthlyMeanTMIN$Month == 
##     "06", ])
## 
## Residuals:
##     Min      1Q  Median      3Q     Max 
## -5.3810 -1.2987  0.0021  1.4953  5.1190 
## 
## Coefficients:
##               Estimate Std. Error t value Pr(>|t|)    
## (Intercept) -23.090329  14.389683  -1.605    0.112    
## YEAR          0.037559   0.007318   5.132 1.47e-06 ***
## ---
## Signif. codes:  0 '***' 0.001 '**' 0.01 '*' 0.05 '.' 0.1 ' ' 1
## 
## Residual standard error: 2.095 on 97 degrees of freedom
## Multiple R-squared:  0.2136, Adjusted R-squared:  0.2055 
## F-statistic: 26.34 on 1 and 97 DF,  p-value: 1.47e-06

## 
## Call:
## lm(formula = TMIN ~ YEAR, data = MonthlyMeanTMIN[MonthlyMeanTMIN$Month == 
##     "07", ])
## 
## Residuals:
##     Min      1Q  Median      3Q     Max 
## -5.7796 -1.6199 -0.1003  1.7751  8.3318 
## 
## Coefficients:
##               Estimate Std. Error t value Pr(>|t|)    
## (Intercept) -14.601690  17.259527  -0.846 0.399631    
## YEAR          0.035276   0.008777   4.019 0.000116 ***
## ---
## Signif. codes:  0 '***' 0.001 '**' 0.01 '*' 0.05 '.' 0.1 ' ' 1
## 
## Residual standard error: 2.516 on 97 degrees of freedom
## Multiple R-squared:  0.1427, Adjusted R-squared:  0.1339 
## F-statistic: 16.15 on 1 and 97 DF,  p-value: 0.0001155

## 
## Call:
## lm(formula = TMIN ~ YEAR, data = MonthlyMeanTMIN[MonthlyMeanTMIN$Month == 
##     "08", ])
## 
## Residuals:
##     Min      1Q  Median      3Q     Max 
## -4.8051 -2.0273 -0.2833  2.0279  7.7257 
## 
## Coefficients:
##               Estimate Std. Error t value Pr(>|t|)    
## (Intercept) -10.458953  17.113818  -0.611 0.542521    
## YEAR          0.033074   0.008702   3.801 0.000251 ***
## ---
## Signif. codes:  0 '***' 0.001 '**' 0.01 '*' 0.05 '.' 0.1 ' ' 1
## 
## Residual standard error: 2.512 on 98 degrees of freedom
## Multiple R-squared:  0.1285, Adjusted R-squared:  0.1196 
## F-statistic: 14.45 on 1 and 98 DF,  p-value: 0.0002506

## 
## Call:
## lm(formula = TMIN ~ YEAR, data = MonthlyMeanTMIN[MonthlyMeanTMIN$Month == 
##     "09", ])
## 
## Residuals:
##     Min      1Q  Median      3Q     Max 
## -5.1896 -1.8215  0.0521  1.6751  9.3185 
## 
## Coefficients:
##              Estimate Std. Error t value Pr(>|t|)    
## (Intercept) -24.61825   18.66327  -1.319     0.19    
## YEAR          0.03916    0.00949   4.127 7.72e-05 ***
## ---
## Signif. codes:  0 '***' 0.001 '**' 0.01 '*' 0.05 '.' 0.1 ' ' 1
## 
## Residual standard error: 2.739 on 98 degrees of freedom
## Multiple R-squared:  0.1481, Adjusted R-squared:  0.1394 
## F-statistic: 17.03 on 1 and 98 DF,  p-value: 7.723e-05

## 
## Call:
## lm(formula = TMIN ~ YEAR, data = MonthlyMeanTMIN[MonthlyMeanTMIN$Month == 
##     "10", ])
## 
## Residuals:
##     Min      1Q  Median      3Q     Max 
## -6.1782 -1.7592 -0.3477  1.8006  6.4333 
## 
## Coefficients:
##              Estimate Std. Error t value Pr(>|t|)   
## (Intercept) -9.124970  18.153567  -0.503  0.61635   
## YEAR         0.028483   0.009232   3.085  0.00265 **
## ---
## Signif. codes:  0 '***' 0.001 '**' 0.01 '*' 0.05 '.' 0.1 ' ' 1
## 
## Residual standard error: 2.647 on 97 degrees of freedom
## Multiple R-squared:  0.08936,    Adjusted R-squared:  0.07997 
## F-statistic: 9.518 on 1 and 97 DF,  p-value: 0.002651

What can explain these contrary trends? It is a phenomenon termed “coastal cooling” that occurs along the entire state of California. When inland temperatures increase due to increased concentration of greenhouse gases, coast to inland pressure and temperature gradients are affected. This causes increased sea breeze frequency, intensity, and/or duration. Therefore, the inland temperature increase from greenhouse gases leads to a reverse reaction of cooling summer maximum temperatures in coastal regions.

Professor Robert Bornstein of San Jose State University, who is at the forefront of research for coastal cooling in California, has conducted research that reflects the temperature changes Ojai is experiencing. Observations done at both regional and global scales have already found that since the mid-1970s, the rate of asymmetric warming has accelerated. Asymmetric warming occurs when TMin increases at a faster rate than TMax, such as what Ojai’s temperature data reveals. Bornstein’s data collection for the South Coast Air Basin area and San Francisco Bay Area show that both TMin and TMax have increased for inland sites, while coastal sites in those regions have simultaneously had rising TMins and falling TMaxes. This discrepancy leads to what appears to be an unchanging average temperature. However, relying solely on the average temperature trendline obscures the descreasing temperatures caused by coastal cooling. This neglects the disruption of the seasonal cycle and potentially harmful impacts. Bornstein’s study concluded that 1970-2005, the overall trend in degrees Celsius per decade for the South Coast Air Basin was: TMin +0.16 for inland areas, +0.28 for coastal areas TMax +0.32 for inland areas, -0.30 for coastal areas

INSERT GRAPHS/PICTURES FROM BORNSTEIN’S PRESENTATION

So what does this all mean? Amidst the warnings of global warming’s destructive impacts on both humans and the environment, any news of cooling should sound like a blessing. The phenomenon of coastal cooling does mean less energy needed for home and building cooling and lower rates of human heat stress.

However, for a town like Ojai rich in organic agriculture, this cooling trend could take a big toll. INFO ABOUT AGRICULTURE - seasons, cool days